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o : 

, The properties of a hard-sphere fluid in contact with hard spherical and cyhndrical waUs are 

fSJ . studied. Rosenfeld's density functional theory (DFT) is applied to determine the density profile 

» I ' and surface tension 7 for wide ranges of radii of the curved walls and densities of the hard-sphere 

CIh, fluid. Particular attention is paid to investigate the curvature dependence and the possible existence 

■ of a contribution to 7 that is proportional to the logarithm of the radius of curvature. Moreover, 

by treating the curved wall as a second component at infinite dilution we provide an analytical 
expression for the surface tension of a hard-sphere fiuid close to arbitrary hard convex walls. The 
agreement between the analytical expression and DFT is good. Our results show no signs for the 
existence of a logarithmic term in the curvature dependence of 7. 
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I. INTRODUCTION 



The properties of inhomogeneous fluids with spherical or cylindrical geometries have been investigated by many 
authors in recent years [1-9] because curvature plays an important role in many physical situations, including the 
^ , following examples. In the context of nucleation phenomena [10-12] one considers the spontaneous formation of 
spherical droplets in a super-saturated vapor. The variation of the local curvature of confinements leads to a more 
complicated form of the depletion forces in colloidal suspensions [13-15]. Wetting or critical adsorption of a fluids 
on curved substrates, such as spherical or cylindrical colloidal particles may lead to the formation of liquid bridges 
if two colloids come close to each other, and subsequently to flocculation [4,7,16-19]. In studies of the hydrophobic 
(~| . effect one is concerned with the cost in free energy to immerse a hydrophobic particle into a liquid [20,21]. Despite 
O " these efforts, the properties of curved fluid interfaces are not understood as fully as the planar one. Whether one 
, O , . deals with a drop of liquid surrounded by its vapor or with a big spherical particle immersed into a solvent of small 
particles, the statistical-mechanical description of such systems generates additional difficulties [22,23] which call for 
I , further studies [24-26]. 

^ • Within this context, here we consider the case of a hard-sphere fluid in the grand canonical ensemble, characterized 
C \ by the radius R of the spheres and the bulk number density p or equivalently the bulk packing fraction r] = (47r/3) p, 
■ in contact with a hard spherical or cylindrical wall. These curved walls can either be viewed as an external potential 
' exerted on the hard-sphere fluid, 

§: -{!!;, :<S, w 

' where the subscript i = s, c denotes the spherical and cylindrical symmetry, respectively, or as the surfaces of a second 
^ . {big) species of radius uf'^ ■ The radii Ri entering the external potential in Eq. (1) are related to the actual radii R^^ 
via i?i = i?f^ -hi? (see Fig. 1). 

^ The grand canonical potential of the system fli can be divided into a bulk part O^""^ = —pVi, where p is the 

Q , pressure of the hard-sphere fluid and Vi is the volume accessible to the system, and an excess (surface) part fi^^: 
O ■ ^f-^^ + ^T- The surface tension 7i is readily defined as the excess grand potential per unit area, 



C)ex Q Cibulk 

><: 



where Af"^ is the surface area of the curved wall. While fli is uniquely defined this division into bulk and surface 
excess parts is not unique. Depending on the choice of the dividing surface the surface tension can even change sign. 
As indicated by the superscript (6) in Eq. (2) we use the actual wall corresponding to the radius Rf"^ as the dividing 
interface. 

Within density functional theory (DFT) , the grand potential Hi of the system can be expressed as 

n,[p] = T[p] + I d\p{v){vr\v) - p) , (3) 
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where J^[p\ ~ T^x [p] + ^id \p\ is the intrinsic free energy functional that can be split into an ideal gas contribution 
Tid [p] and an excess (over the ideal gas) term T^x [p] • The only explicit dependence of [p] on the external potential 
originates from the second term in Eq. (3). The equilibrium density profile p(r) satisfies the Euler-Lagrange equation 

(5p(r) 

For the equilibrium density profile p(r) the density functional reduces to the grand potential of the system, i.e., 
fij = f2j[p(r)], and we can rewrite Eq. (2) as 

n,[p(r)]+py, 

7i = • (5) 

From Eqs. (1), (3), and (4) it follows that the change in the grand potential due to an infinitesimal change of the 
radius i? •''^ of the curved wall, at constant chemical potential /x and temperature T, is given by 

The first term on the right-hand side of Eq. (6) vanishes by virtue of Eq. (4) while the second term, for the external 
potential given in Eq. (1), is equal to l3~^Ai p{Rf), where piRf) denotes the contact density of the hard-sphere fluid 
at the curved wall. This leads to the sum rules [24] 

/ d^i \ _ J ATrR'^p{R+), for a spherical wall , . 

\dR^ ) j,~ \ 27ri?cip(-R+), for a cylindrical wall. ^ ' 

In the case of a cylindrical wall we consider the thermodynamic limit lim {Q.c/L) so that effects due to a finite length 

L— »oo 

L of the cylinder drop out. Wc note that the smu rules in Eq. (7) arc valid for all one-component fluids in contact with 
hard spherical or cylindrical walls and are satisfied by all density functionals within weighted-density approximation 
(WDA) [2]. 

Using the actual radius R^^ also as the radius of the dividing interface (see Fig. 1), the grand potential = 
—pVi + ^iAf'' of the system is separated into a bulk and surface term with the accessible volume 

Vi = Vi{R)-Vi{Rf^), (8) 

where Vs{r) — , K(f') = Trr'^L, and i? is a macroscopicly large radius considered in the thermodynamical limit, 
and the surface areas 



^(b) _ I 47r(i?i''^)^, for a spherical wall 
* I 2-kR^^L, for a cylindrical wall . 



(9) 



The sum rule in Eq. (7) can then be expressed as 



for a spherical wall and 



piRc)-m=Pp+ — -n^+P{—-n^] (11) 



itc ric \ Ottc 
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for a cylindrical wall. Equations (10) and (11) indicate that different versions of DFTs might lead to different contact 
values near curved walls even if the underlying equation of state is the same, which leads to identical contact values 
at a hard planar wall. As pointed out by J.R. Henderson [26], these sum rules have interesting implications. If p{R'l) 

is analytic as a function of Rr^ = (i?f ^ + R^^ one has 
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with Co = I3p for both the spherical and the cylindrical wall and expansion coefficients cf ^ for j = 1,2, . . .. Note 
that the expansion in Eq. (12) is assumed to converge for arbitrary non-zero values of iij. This assumption will be 
corroborated by our numerical results in Sec. III. With Eq. (12) the solution of the differential equation (10) can be 
written as 



(1) 



1 + 
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(2) 
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(13) 



(1) 



cP + 2ppR and cP 



where Ds is an integration constant, C. 
solution of the differential equation (11) can be written as 

4'^ ln(4''V^) 



cl^^ +ci^^i? + /Jpi?2. Likewise, with Eq. (12) the 
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where D, 



is an integration constant and Ci^' = ci^^ 



f3pR. In the case of the spherical wall 7^ {f6s'' ) has the form 



/374^«?^)=/?7a(oo) 
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where 7^ (00) = /3 ^ %- + pR is the surface tension of the planar wall and 



cf^ +cPr + PpR'' 
cP + 2(3pR 



(16) 



is the so-called Tolman length [27], which plays an important role in nucleation theory [4,10-12]. We note that the 

value of the Tolman length depends on the choice of the dividing surface; here the dividing surface has the radius 

R^s'' ■ If these logarithmic terms in the expansion of the surface tensions do not vanish, the concept of the Tolman 
length apparently breaks down for cylindrically curved and generally shaped walls [25,26]. 

The analyticity of 7(-R) in terms of is of considerable importance because it provides the basis for the so-called 
Helfrich expansion of the surface free energy of arbitrarily curved surfaces in powers of the principal curvatures [28,29]. 
Although this approach appears to be very useful for describing membranes [30] it has been shown recently that the 
presence of long-ranged dispersion forces destroys the aforementioned analyticity for fluid interfaces [4,31,32] and 
prevents the use of a Helfrich-typc interface Hamiltonian for such systems. The possible occurrence of logarithmic 
terms in Eqs. (13) and (14) following from the simple expansion (12) raises the issue whether the absence of analyticity 
is not only due the presence of long-ranged forces but is an intrinsically geometrical effect which shows up even for 
surface free energies of hard-sphere fluids. The present study aims at clarifying this conceptual point. 

To this end we investigate the curvature dependence of the surface tension as well as relationships between ther- 
modynamic and local properties of the hard-sphere fluids following two different routes. First, in Sec. II the surface 
tension is determined from a bulk theory in which the curved wall is considered as the surface of an additional com- 
ponent at infinite dilution. Second, a numerical approach based on the minimization of the Rosenfeld fundamental 
measure theory (FMT) free energy functional [33,34] is presented in Sec. III. In this approach, the hard-sphere fluid 
is considered to be exposed to the external potential of Eq. (1). 



II. BULK THEORY 



We start by considering the bulk of a one component hard-sphere fluid in a volume Vtot with Vtot — in the 
thermodynamic limit. The grand potential of this system is given by flo = —pVtot, where p is the pressure of the 
hard-sphere fluid. A single particle of a second component b is inserted into this system. In view of the following 
discussion the only restriction for this particle is that it must be hard and convex. The grand potential of the new 
system is 



(17) 



where Arj measures the change of the grand potential due to the insertion of a single particle. Ail equals the one- 
particle direct correlation function [14,15,35] which in the bulk limit considered here is simply the excess chemical 
potential of species b in the dilute limit pb — > 0, i.e., 

f3An = -4^^ . (18) 

Although Eqs. (17) and (18) are formally exact it is important to notice that this approach within a bulk theory is 
expected to only work reliably if the considered fluid is in a single phase. In the case of phase separation, which could 
lead to wetting or drying of the big sphere, this approach most likely would fail [36]. 

We evaluate Eq. (18) within the framework of Rosenfeld's Fundamental Measure Theory (FMT). Within this 
approach [33,34] the excess free energy functional (over the ideal gas) is given by 

PJ^ex = J d^r $({n„}) , (19) 

where Ua. denote weighted densities 

N 

n« (r) = ^ / d'r'pi (r') w^'^ (r - r') , (20) 

with geometrical weight functions Wa of species j = 1, . . . ,N. In Rosenfeld's approach there are six different weighted 
densities, four are scalar and two are vector-like. For a general hard body j one has [37,38] 

«;(j')(r) = e(|r-R(^)|) , 
«;^^'(r) =5(|r-R(^')|) , 
w(^](r) =n(^') 5(|r-R(^')|) , 

U)i- ^ ^^^^ U)f ^ 

wl^](r) = ^w(^](r) , (21) 

where R^^^ — R^J\6,ip) is the radius vector to the surface and n(j') is its outward unit vector normal to the surface 
at point R*^-'^ of particle j. H^^^ and K^^^ are the integrated mean curvature and integrated Gaussian curvature of 
the particle j, respectively [37,38]. 

In the bulk limit for which the density profiles pj{r) are constant, as considered in this section, the vector weighted 
densities vanish, and Eq. (18) takes the form 

-4^ = T— Cs + j^ C2 + T— Ci + T— Co, 22 
ans on2 on\ ouq 

with characteristic functions Ci, « = 0, . . . , 3 of the shape of the particle of species b. 

Assuming the general functional of the given by Eqs. (19)-(21), one has for a general convex hard particle [37,38] 

C3 = V^''\ (23) 
C2 = AW, (24) 

Cr=AW^, (25) 

Co=AW^, (26) 

where the shape of the particle is specified by its volume its surface area A'^''\ and integrated curvatures i?^''^ 

and K^^l 



and 
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Within FMT there are several expressions for the excess free energy density For the given problem we have 
chosen the original Rosenfeld form [33] 



nin2 - 111 • n2 - 3n2n2 • 112 

1 -ng ^ 247r(l - n^^ 



$({na}) = -no ln(l - ns) + ^ , J + ' ^ : • (27) 



The partial derivatives of the excess free energy density with respect to the weighted densities in the dilute limit of 
species b depend only on the packing fraction r] and the radius R of the hard-sphere fluid: 

a$ 377 l + ri + Tf_ 



dn2 87ri?2(l-?7)2 

a$ 3?? 



am R{1 - T]) 

and 



, (30) 



= -ln{l-r,)^<l>o. (31) 

ono 

In the above equations ppy and are the Percus-Yevick compressibility equation of state and the scaled-particle 

theory (SPT) expression for the surface tension of a hard-sphere fluid at a planar hard wall [39], respectively. 

With the characteristic functions in Eqs. (23)-(26) and the partial derivatives given by Eqs. (28)-(31) one flnds 
within this approach for the grand potential of a hard-sphere fluid surrounding any convex hard particle b 

/ TT(b) T^'(b) \ 

= -PppviVtot - V^'y) + At") i^p^l^PJ + ^ ^1 + ^ 0oj . (32) 

By choosing the surface of the convex particle as the dividing interface, the bulk contribution to the grand potential 
is simply flo = —p^^{Vtot — V^''^) and the surface tension follows directly as 

/37(//(^\i^W) = p^SPT + ^L_^^ + !L_^^, (33) 

The particular choice of the dividing interface underlying Eq. (33) is motivated by the simplicity of the resulting 
expression for the surface tension. Other choices of the dividing interface would lead to significantly more complex 

expressions and the elegant, simple form of Eq. (33) would be lost. Obviously within this bulk approach it is impossible 
to pick up a contribution to the surface tension that is proportional to the logarithm of the radius of curvature. 

The surface tension in Eq. (33) is always positive for a non-zero packing fraction 77 and its change for small curvature 
due to an increase in curvature is described by the Tolman length St which within this approach is defined by 

7(//W, /f(^)) = ^'h^ (1 - 2H^»^6t + O(i^W)) , (34) 

hape of the convex particle with 



so that it is independent of the shape of the convex particle with 
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Thus within this approach we find that the generalized Tolman length for a general hard convex cavity is negative, 
provided the surface of this hard cavity is chosen as the dividing surface. The grand potential (Eq. (32)) and the 
surface tension (Eq. (33)) have been derived for a general convex hard wall and are related to the Helfrich Hamiltonian 
[28] . In order to be able to calculate the contact density of a hard-sphere fluid on the basis of the sum rule in Eq. (7) 

one has to evaluate the expression in Eq. (33) for special cases. For a big spherical particle the mean curvature 

ni''^ = l/R^s^ and the Gaussian curvature K^^ = l/{Ri''^)'^. By combining Eqs. (33) and (10) within the present 
bulk theory one finds for the contact value of the density profile of a hard-sphere fluid at a big hard sphere: 
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For finite curvatures the contact density is smaller than in the planar limit, which in the present approach is given 
by the Percus-Yevick compressibility pressure f3ppY- The difference between the contact density at finite curvature 
and its planar wall limit is a quadratic function in 1/Rs- 

For a cylindrical wall we have to proceed more carefully. Since Eq. (33) was derived assuming a finite convex 
body we take the limit L — > oo only after the cylinder is inserted into the bulk system. In this limit the mean 
curvature Hc''^ = l/{2Ri!'^) and the Gaussian curvature Kc''^ = 0. The contact density of a hard-sphere fluid close to 
a cylindrical wall follows from Eq. (33) and the sum rule in Eq. (11) as 

R 9r;2 / (1 + r?) E \ 

which is linear in l/i?c- For large radii of curvature Ri oc the contact densities at a big sphere p{Rf) and at a 
cylinder p{R+) are related by 2p{R+) - p{R+) = Pppy + 0{^). 



III. NUMERICAL RESULTS 

As mentioned in the Introduction an alternative approach for calculating the surface tension and the contact 
density of a hard-sphere fluid at a spherical or cylindrical wall consists of minimizing the density functional, Eq. (3), 
numerically. Within this approach the curved wall enters the problem via the external potential specifled in Eq. (1). 

The solution of the corresponding Eulcr-Lagrange equation yields the inhomogcneous equilibrium density profile p{r). 

Prom p(r) both the contact densities p{R^) = p{r = R^) and the surface tensions ji{Rf ^), Eq. (5), follow directly. 
In order to verify that the numerical DFT results for the contact density and the surface tension (from minimization) 
arc consistent and to probe the numerical accuracy of our DFT calculation, we confirmed numerically that the sum 
rule in Eq. (7) is satisfied very accurately (within the range of four to six significant digits depending on the packing 
fraction). 

Although the minimization of the density functional and the bulk route discussed in Sec. II are in principal equivalent 
approaches, they differ in practice because the actual density functional for the hard-sphere fluid is an approximation. 
Within an approximate treatment it is expected that the numerical solution of the Euler-Lagrange equation, which 
is based on a variational principle, is more accurate than the aforementioned bulk approach. Nevertheless, for the 
hard-sphere fluid which, apart from freezing, exhibits only a single phase both routes are expected to be in qualitative 
or even semi-quantitative agreement with each other. 

First we present in Fig. 2 our numerical results for the surface tension at a spherical wall as obtained by full 
minimization of the Rosenfeld functional for values of ry in the range between ?7 = 0.1 and r] = 0.42 and with radii 
of curvature Ri^^ in the range Ri'''' = R and = lOOOOi?. The DFT results arc denoted by symbols and the 

predictions of the bulk theory by full lines. For small values of rj we flnd both results to be in excellent, quantitative 
agreement for all curvatures, see Fig. 2(a). Close to the planar wall limit the level of quantitative agreement reduces 
signiflcantly for larger packing fractions, while it remains good for large curvatures. This deviation between the two 
different routes is expected because in the limit Ri''^ — » oo the bulk route reduces to the SPT prediction of the surface 
tension at a planar wall, which is known to overestimate the surface tension of a hard-sphere fluid at a planar hard 
wall [40]. Since SPT obtains the planar wall limit from extrapolating the behavior at small cavity sizes [39] this 
deviation is no surprise. Recently it was suggested that SPT can be re-optimized by taking the planar wall limit into 
account more accurately [36]. 

Even for packing fractions at which the level of quantitative agreement is moderate (see Fig. 2(b)), there is a good 
qualitative agreement between the DFT results and our bulk approach. This implies that the curvature dependence 
of our numerical data of 7s can be best approximated by a polynomial quadratic in l/R^'' with no evidence for a 
logarithmic term. 

The DFT results for the surface tension of a hard-sphere fluid in contact with a hard cylindrical wall are shown in 
Fig. 3. Again the quantitative agreement between the DFT results (symbols) and predictions of the bulk theory (full 
lines) are excellent for low packing fractions and fair at higher values of ry. The curvature dependence of the DFT 

data is captured very well by a function linear in 1/ R^c^ as predicted by the bulk route. 

A remarkable feature of Eq. (33) is that the ratio of the coefficients corresponding to the term of the leading order 

in 1/Ri for spherical and cylindrical symmetry is exactly 2 (since H^J''^ = 2Hlt^). This behavior can also be found in 
the numerical results of DFT calculations, i.e., based on the inhomogcneous density profiles p{r). In Table I wc show 

values of least-square fits assuming that f3jsiRs'^) = ^ + -jej + "l,)^ and jdRc^^) = Oc*''' + -jbjj i-e., the simplest 
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formulae to fit our DFT data. To a good approximation the ratio ai /oc equals 2 for all densities. This ratio agrees 
with the predictions of the bulk theory; a^^^ = a^c^ is the planar wall surface tension. The Tolman length (compare 

Eq. (34)) is given by Si^^ = — ai^'/(2ai°^) and — —a^c^/a^c^. In Fig. 4 we plot (5^*' (diamonds) and (5^^ (circles) 
together with the prediction of the bulk theory Eq. (35) (full line) as a function of the packing fraction r). 

Thus neither the bulk theory nor the numerical DFT results provide any hints for the occurrence of logarithmic 
terms in the expression of the surface tension of hard-sphere fluids in contact with curved hard surfaces. 

For the contact density one knows a priori that for all values of rj both the DFT route and the bulk approach yield 
the same contact value p{Rf oo) = (3p of the density profile of a hard-sphere fluid in the limit of the planar wall. 
For finite curvature we find that the contact densities obtained numerically from DFT (symbols) and analytically 
from our bulk theory for 7(i?i''^) (full lines) are in very good quantitative agreement for all values of rj, as shown in 
Fig. 5 for spherical walls and in Fig. 6 for cylindrical walls. Only for very large curvatures the actual DFT contact 
densities are slightly overestimated by the bulk theory. 

For the spherical wall in Fig. 5 we also plot the contact density obtained from an empirical fit to simulation data 
(dotted lines) by Degreve et al. [6]. In terms of Rj^ this fit interpolates linearly between the planar wall limit 
within the Carnahan-Starling theory and the contact value of the pair correlation function g{r) corresponding to 
R/Rs = 0.5. At high packing fractions this fit deviates from our results in the planar wall limit due to the difference 
between the Carnahan-Starling equation of state and the Percus-Yevick compressibility pressure which underlies both 
our approaches. At very high packing fractions there are additional deviations because the simple linear fit does not 
capture the actual higher order terms in R~^. Nonetheless, the overall agreement between our results and this fit is 
good. 

IV. CONCLUSIONS AND SUMMARY 

Our analysis of the local structure and the thermodynamics of the hard-sphere fluid near hard curved substrates 
has lead to the following main results: 

1. We have studied the curvature dependence of the surface tension of a hard-sphere fluid at hard spherical and 
cylindrical walls (Fig. 1) obtained by minimizing the Rosenfeld fundamental measure density functional. We have 
found no indications of logarithmic singularities in the expansion of the surface tension in terms of curvatures 
of spherical or cylindrical walls. The surface tension appears to be an analytic function of curvature. 

2. Based on the Rosenfeld fundamental measure theory we have derived an analytical expression for the surface 
tension of a hard-sphere fluid close to hard, arbitrarily curved, convex walls in terms of their integrated mean and 
Gaussian curvatures, H^''^ and K^''^ (see Eq. (33)). This approach also does not render logarithmic singularities 
of the surface tension as a function of the radii of curvature. There is good agreement between the results of 
the density functional theory and the bulk theory which treats the curved wall as the siirface of a particle of a 
second component in infinite dilution (Figs. 2 and 3). The Tolman length as a function of the packing fraction 
of the hard-sphere fluid is shown in Fig. 4. 

3. Based on the sum rule for curved substrates we have obtained expressions for the contact values of the density 
profile of a hard-sphere fluid close to spherical (Eq. (36)) and cylindrical (Eq.(37)) walls. A comparison with 
the density functional results indicates, that the analytical expressions of the contact values are very reliable; 
only at rather high packing fractions they slightly overestimate the density functional values (Figs. 5 and 6) for 
large curvatures. 

4. The ratio of the leading order terms in the expansion of the surface tension as a function of curvature, obtained 
from full minimization of the fundamental measure functional, for cylindrical and spherical symmetry equals 2 
(see Table I). This is in agreement with the aforementioned bulk theory (Eq. (33)). This also agrees with the 
general feature of the Helfrich theory [28], in which for a sphere the contribution to the surface tension linear 

in the curvature is twice as large as for a cylinder. 

In Ref. [41] the arguments in favor of and against the existence of logarithmic singularities in the curvature depen- 
dence of the surface tension of fluids have been reviewed. Our present results indicate that for a hard-sphere fluid 
near hard curved substrates such logarithmic singularities either do not exist or their influence is negligible. 

It would be of interest to investigate the issue of analyticity of the surface tension via a diagrammatic expansion of 
the contac;t value of the density profile [42,43]. Within this approach our preliminary calculations show, that in the 
lowest order in density the diagrammatic expansion of the contact value of the density profile indeed leads to a surface 
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tension analytical in curvature. However, the extension of this analysis to higher orders in densities is technically 
involved and is postponed to future work. 

Thus our findings indicate that non-analyticities of the curvature dependence of surface tensions arise only via 
dispersion forces acting on the fluid particles [4] or via the onset of drying or wetting transitions on curved substrates 

[44] 
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FIG. 1. Geometry of a hard-sphere fluid in contact with a 
Rf =Ri-R and volume V^^\ While the density profile p(r) 
surface with radius Rf"^ as dividing surface for the calculation 



hard spherical {i = s) or cylindrical {i = c) wall with radius 
attains its contact value at r = iij, we choose the wall, i.e., the 
of the surface tension. 
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FIG. 2. The surface tension of a fluid of hard spheres with radius R near a hard spherical wall of radius Rl ■ Symbols denote 

results obtained from direct minimization of the Rosenfeld functional, whereas the full lines are predictions of the bulk theory 
(see Eq. (33)). Small packing fractions (77 = 0.1, . . . , 0.25) are shown in (a) and large packing fractions (77 = 0.3, . . . , 0.42) in (b). 
The agreement between both routes is excellent for low packing fractions. Towards high values of r] the level of quantitative 
agreement decreases, while the bulk theory still reproduces qualitatively the DFT results. 
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FIG. 3. The surface tension of a fluid of hard spheres with radius R near a hard cyUndrical wall of radius Rc ■ Symbols denote 

results obtained from direct minimization of the Rosenfeld functional, whereas the full lines are predictions of the bulk theory 
(see Eq. (33)). Small packing fractions (77 = 0.1, . . . , 0.25) are shown in (a) and large packing fractions (77 = 0.3, . . . , 0.42) in (b). 
The agreement between both routes is excellent for low packing fractions. Towards high values of r] the level of quantitative 
agreement decreases, while the bulk theory still reproduces qualitatively the DFT results. 
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FIG. 4. The Tolman length 5t of a fluid of hard spheres with radius R at curved hard walls as a function of the packing 
fraction rj. Diamonds and circles denote results for spherical and cylindrical walls, respectively, while the full line correspond 
to the prediction of the bulk theory (Eq. (35)). 
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FIG. 5. The contact values of the density profile of a fluid of hard spheres with radius _R at a hard spherical wall with radius 
— Rs — R sls obtained from direct minimization of the Rosenfeld functional (symbols), from the bulk theory (Eq. (36)) 
(solid lines), and from a semi-empirical parameterization [6] (dotted lines). The results for the packing fraction rj = 0.2 are 
given in (a) and for r] = 0.4 and r] = 0.42 in (b). Even for high packing fractions the agreement between the DFT results and 
the predictions of the bulk theory is very good. 
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R/R^ 

FIG. 6. The contact values of the density profile of a fluid of hard spheres with radius J? at a hard cylindrical wall with 
radius R^^ = Rc — R as obtained from direct minimization of the Rosenfeld functional (symbols) and from the bulk theory 
(Eq. (37)) (solid lines). The results for the packing fractions r; = 0.1 and rj = 0.2 arc given in (a) and for q — 0.4 and r; — 0.42 
in (b). Even for high packing fractions the agreement between the DFT results and the predictions of the bulk theory is very 
good. 
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TABLE I. Values of a^^\ o.i^\ ai"^' and ai^^ as obtained from least square fit to DFT data of the surface tension of a 
hard-sphere fluid in contact with a spherical and cylindrical hard wall for various values of the packing fraction r]. Note that 

Eq. (33) predicts 4iy = 2 for all r). 
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